A scheme for radiation pressure and photon diffusion with the M1 closure in RAMSES-RT

We describe and test an updated version of radiation-hydrodynamics (RHD) in the RAMSES code, that includes three new features: i) radiation pressure on gas, ii) accurate treatment of radiation diffusion in an unresolved optically thick medium, and iii) relativistic corrections that account for Doppler effects and work done by the radiation to first order in v/c. We validate the implementation in a series of tests, which include a morphological assessment of the M1 closure for the Eddington tensor in an astronomically relevant setting, dust absorption in a optically semi-thick medium, direct pressure on gas from ionising radiation, convergence of our radiation diffusion scheme towards resolved optical depths, correct diffusion of a radiation flash and a constant luminosity radiation, and finally, an experiment from Davis et al. of the competition between gravity and radiation pressure in a dusty atmosphere, and the formation of radiative Rayleigh-Taylor instabilities. With the new features, RAMSES-RT can be used for state-of-the-art simulations of radiation feedback from first principles, on galactic and cosmological scales, including not only direct radiation pressure from ionising photons, but also indirect pressure via dust from multi-scattered IR photons reprocessed from higher-energy radiation, both in the optically thin and thick limits.Comment: 25 pages, 13 figures, accepted for publication in MNRAS. Revised to match published versio

Globular Cluster Formation in the Virgo Cluster

Metal poor globular clusters (MPGCs) are a unique probe of the early universe, in particular the reionization era. Systems of globular clusters in galaxy clusters are particularly interesting as it is in the progenitors of galaxy clusters that the earliest reionizing sources first formed. Although the exact physical origin of globular clusters is still debated, it is generally admitted that globular clusters form in early, rare dark matter peaks (Moore et al. 2006; Boley et al. 2009). We provide a fully numerical analysis of the Virgo cluster globular cluster system by identifying the present day globular cluster system with exactly such early, rare dark matter peaks. A popular hypothesis is that that the observed truncation of blue metal poor globular cluster formation is due to reionization (Spitler et al. 2012; Boley et al. 2009; Brodie & Strader 2006); adopting this view, constraining the formation epoch of MPGCs provides a complementary constraint on the epoch of reionization. By analyzing both the line of sight velocity dispersion and the surface density distribution of the present day distribution we are able to constrain the redshift and mass of the dark matter peaks. We find and quantify a dependence on the chosen line of sight of these quantities, whose strength varies with redshift, and coupled with star formation efficiency arguments find a best fitting formation mass and redshift of $\simeq 5 \times 10^8 \rm{M}_\odot$ and $z\simeq 9$. We predict $\simeq 300$ intracluster MPGCs in the Virgo cluster. Our results confirm the techniques pioneered by Moore et al. (2006) when applied to the the Virgo cluster and extend and refine the analytic results of Spitler et al. (2012) numerically.Comment: 13 Pages, 13 Figures, submitted to MNRA

Infall near clusters of galaxies: comparing gas and dark matter velocity profiles

We consider the dynamics in and near galaxy clusters. Gas, dark matter and galaxies are presently falling into the clusters between approximately 1 and 5 virial radii. At very large distances, beyond 10 virial radii, all matter is following the Hubble flow, and inside the virial radius the matter particles have on average zero radial velocity. The cosmological parameters are imprinted on the infall profile of the gas, however, no method exists, which allows a measurement of it. We consider the results of two cosmological simulations (using the numerical codes RAMSES and Gadget) and find that the gas and dark matter radial velocities are very similar. We derive the relevant dynamical equations, in particular the generalized hydrostatic equilibrium equation, including both the expansion of the Universe and the cosmological background. This generalized gas equation is the main new contribution of this paper. We combine these generalized equations with the results of the numerical simulations to estimate the contribution to the measured cluster masses from the radial velocity: inside the virial radius it is negligible, and inside two virial radii the effect is below 40%, in agreement the earlier analyses for DM. We point out how the infall velocity in principle may be observable, by measuring the gas properties to distance of about two virial radii, however, this is practically not possible today.Comment: 7 pages, 3 figures, to appear in MNRA

Modelling CO emission from hydrodynamic simulations of nearby spirals, starbursting mergers, and high-redshift galaxies

We model the intensity of emission lines from the CO molecule, based on hydrodynamic simulations of spirals, mergers, and high-redshift galaxies with very high resolutions (3pc and 10^3 Msun) and detailed models for the phase-space structure of the interstellar gas including shock heating, stellar feedback processes and galactic winds. The simulations are analyzed with a Large Velocity Gradient (LVG) model to compute the local emission in various molecular lines in each resolution element, radiation transfer and opacity effects, and the intensity emerging from galaxies, to generate synthetic spectra for various transitions of the CO molecule. This model reproduces the known properties of CO spectra and CO-to-H2 conversion factors in nearby spirals and starbursting major mergers. The high excitation of CO lines in mergers is dominated by an excess of high-density gas, and the high turbulent velocities and compression that create this dense gas excess result in broad linewidths and low CO intensity-to-H2 mass ratios. When applied to high-redshift gas-rich disks galaxies, the same model predicts that their CO-to-H2 conversion factor is almost as high as in nearby spirals, and much higher than in starbursting mergers. High-redshift disk galaxies contain giant star-forming clumps that host a high-excitation component associated to gas warmed by the spatially-concentrated stellar feedback sources, although CO(1-0) to CO(3-2) emission is overall dominated by low-excitation gas around the densest clumps. These results overall highlight a strong dependence of CO excitation and the CO-to-H2 conversion factor on galaxy type, even at similar star formation rates or densities. The underlying processes are driven by the interstellar medium structure and turbulence and its response to stellar feedback, which depend on global galaxy structure and in turn impact the CO emission properties.Comment: A&A in pres

A scheme for radiation pressure and photon diffusion with the M1 closure in ramses-rt

We describe and test an updated version of radiation-hydrodynamics in the ramses code, that includes three new features: (i) radiation pressure on gas, (ii) accurate treatment of radiation diffusion in an unresolved optically thick medium, and (iii) relativistic corrections that account for Doppler effects and work done by the radiation to first order in v/c. We validate the implementation in a series of tests, which include a morphological assessment of the M1 closure for the Eddington tensor in an astronomically relevant setting, dust absorption in an optically semithick medium, direct pressure on gas from ionizing radiation, convergence of our radiation diffusion scheme towards resolved optical depths, correct diffusion of a radiation flash and a constant luminosity radiation, and finally, an experiment from Davis etal. of the competition between gravity and radiation pressure in a dusty atmosphere, and the formation of radiative Rayleigh-Taylor instabilities. With the new features, ramses-rt can be used for state-of-the-art simulations of radiation feedback from first principles, on galactic and cosmological scales, including not only direct radiation pressure from ionizing photons, but also indirect pressure via dust from multiscattered IR photons reprocessed from higher-energy radiation, both in the optically thin and thick limit

Collapse, outflows and fragmentation of massive, turbulent and magnetized prestellar barotropic cores

Stars and more particularly massive stars, have a drastic impact on galaxy evolution. Yet the conditions in which they form and collapse are still not fully understood. In particular, the influence of the magnetic field on the collapse of massive clumps is relatively unexplored, it is thus of great relevance in the context of the formation of massive stars to investigate its impact. We perform high resolution, MHD simulations of the collapse of hundred solar masses, turbulent and magnetized clouds, using the adaptive mesh refinement code RAMSES. We compute various quantities such as mass distribution, magnetic field and angular momentum within the collapsing core and study the episodic outflows and the fragmentation that occurs during the collapse. The magnetic field has a drastic impact on the cloud evolution. We find that magnetic braking is able to substantially reduce the angular momentum in the inner part of the collapsing cloud. Fast and episodic outflows are being launched with typical velocities of the order of 3-5 km s$^{-1}$ although the highest velocities can be as high as 30-40 km s$^{-1}$. The fragmentation in several objects, is reduced in substantially magnetized clouds with respect to hydrodynamical ones by a factor of the order of 1.5-2. We conclude that magnetic fields have a significant impact on the evolution of massive clumps. In combination with radiation, magnetic fields largely determine the outcome of massive core collapse. We stress that numerical convergence of MHD collapse is a challenging issue. In particular, numerical diffusion appears to be important at high density therefore possibly leading to an over-estimation of the number of fragments.Comment: accepted for publication in A&